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TABLES  FOR  DETERMINATION  OF  FLOW  VARIABLE  GRADIENTS 
BEHIND  CURVED  SHOCK  WAVES 

ABSTRACT 

The  shock  wave  relations  and  the  equations  of  isoenergetic 
two-dimensional  or  axisymmetric  flows  can  be  combined  to  yield 
expressions  for  various  useful  quantities  behind  a  shock  wave  if 
the  curvature  is  known  in  addition  to  the  position  and  slope.  Tables 
are  presented  here  of  computations,  for  7  =  1.4,  of  seven  functions 
of  Mach  number  and  shock  wave  slope,  from  which  the  following  are 
easily  computed:  l)  streamline  curvature,  2)  velocity  gradient  along 
streamline,  3)  angle  between  streamline  and  density  contour,  and  4) 
angle  between  streamline  and  Mach  number  contour.  In  addition,  the 
Mach  number  and  angle  of  inclination  of  streamline  are  listed.  A 
Mach  number  range  from  1.1  to  10  is  covered  by  these  computations. 

A  derivation  of  the  shock  gradient  functions  is  presented;  and 
several  applications  of  the  calculations  are  given,  including  a  set 
of  tables  for  determining  the  slope  of  the  sonic  line. 
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-  rectangular  coordinates 

-  natural  coordinates  normal  and  tangential  respectively 
to  streamline 

l  2  _  2,2  qn1/2 

=  (cos  p  +  g  sin  p)  ' 

-  speed  of  sound  *=  (dp /dp)"^ 

-  shock  gradient  functions 

-  inverse  of  density  ratio  across  shock  =  p^/p 

-  enthalpy 

-  curvature  of  shock  wave 

-  curvature  of  streamline 

-  inverse  of  pressure  ratio  across  shock  =  p^/p 

-  Mach  number 
=  sin  P 

-  pressure 

-  velocity 

-  ratio  of  stagnation  densities  across  shock 

-  entropy 

-  temperature 

-  angle  between  shock  wave  and  x-axis 

-  ratio  of  specific  heats  (=  1.4  for  air) 

-  0  for  2 -dimensional  flow 
1  for  axisymmetric  flow 

-  angle  between  streamline  and  line  of  constant  density 

-  angle  between  streamline  and  x-axis 

-  angle  between  streamline  and  line  of  constant 
Mach  number 
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p  -  density 

a  -  axe  length  along  shock  wave 

Subscripts 

1  -  free  stream 

t  -  stagnation 
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I.  INTRODUCTION 


The  problem  of  determining  the  gradients  of  the  flow  variables 

behind  a  shock  wave  in  isoenergetlc  two-dimensional  flow  has  been 

1  ?  B 

considered  by  several  investigators  ’  .  These  gradients  are  very 

useful;  for  instance,  the  slope  of  the  streamline  at  shock  polar 
points  in  the  holograph  plane  can  be  obtained  from  them  (the  Busemarm 

If 

"hedgehog").  A  particularly  simple  derivation  was  Indicated  by  Sternberg  , 

using  natural  coordinates.  It  is  found  that  the  flow  variable  gradients 

are  proportional  to  the  shock  wave  curvature .  Tables  of  coefficients  for 

2 

streamline  curvature  have  been  computed  by  Thomas  for  a  limited  range  of 
Mach  number. 

The  same  gradients  can  be  computed  for  axi symmetric  flow.  However, 
they  are  now  linear  combinations  of  the  two  shock  curvatures,  K  in  the 

D 

meridianal  plane  and  l/r  in  the  azimuthal  plane.  (The  two-dimensional 
result  in  obtained  simply  by  neglecting  the  l/r  term.)  Because  of  this 
the  general  usefulness  of  these  gradient  functions  is  restricted  in  that 
one  can  no  longer,  as  in  two-dimensional  flow,  obtain  certain  quantities 
(e.g.,  slope  of  sonic  line)  without  having  to  specify  the  values  of  r  and 
K  .  Expressions  for  the  slope  of  these  contours  (density,  pressure,  Mach 
number)  have  been  given  by  Wood  and  Gooderum  . 


(a) 


(b) 


x 


FIGURE  l.l 
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Aside  from  the  results  of  Thomas  cited  above  no  calculations 
of  the  above-mentioned  gradient  quantities  have  to  the  authors’  knowledge 
been  published.  This  report  supplies  such  a  set  of  calculations  ("shock 
gradient  functions")  which  can  be  used  to  determine  significant  flow  in¬ 
formation  from  physical  measurements .  It  is  felt  that  these  computations 
can  also  be  helpful  in  obtaining  a  clearer  general  picture  of  flow  behind 
a  curved  shock  wave.  These  functions  will  find  further  application  in  the 
study  of  flow  of  a  relaxing  gas  behind  shock  waves  at  large  Mach  numbers. 

The  particular  quantities  to  be  obtained  here  are  (using  the  notation 
of  Fig.  l.l):  l)  dq/ds,  the  velocity  gradient  along  the  streamline)  2) 
d©/8s  (=  K^),  the  curvature  of  the  streamline;  5)  ©  +  £,  the  slope  angle 
of  the  isopycnal,  or  density  contour;  and  k)  0  +  1,  the  slope  angle  of  the 
Mach  number  contour.  It  is  assumed  that  the  gas  is  ideal,  and  that  the 
flow  is  non-viscous  and  isentropic  along  streamlines. 


The  above  flow  quantities  are  obtainable  from  expressions  having  the 
following  simplified  forms: 


d©/ds 

(l/^dq/ds 


tan  £ 


-  ^(7,  Mj,  P)  Kg  +  e  F2(7,  Mx,  p)  (l/r) 

=  F5(7,  Mj,  P)  Ks  +  e  Fu(7,  M^  p)  (l/r) 

"CT  f  XI  q)  V  ■  TP  ( TUT  Q  ^  ( ~\  t  ■*+\ 

=  F'5(7,  Mx,  p;  Ks  +  e  F6(7,  Mx,  p)  (l/r) 


(1.1) 


F,(7,  H-l,  P)  Kg  +  e  F J7,  M^  p)  (l/r) 
tan  1  =  F^(7,  ^7  ^)  Ks  +  e  F6(7,  M^  p)  (l/rj 

where  e  is  equal  to  zero  for  two-dimensional  flow  and  one  for  axi symmetric 
flow.  M1  and  q1  are  the  Mach  number  and  velocity,  respectively,  in  the 
free  stream;  tan  P  is  the  slope  of  the  shock  wave;  7  is  the  ratio  of  specific 
heats.  The  angle  0=0  (M^,  p)  is  known  from  the  shock  wave  relations 
(see,  e.g..  Ref.  6), 
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The  gradient  functions  F^,  F2,  . are  derived  in  detail  in 
Sectionll.  These  functions  for  the  most  part  are  long  and  complicated 
expressions  in  and  p  which  do  not  lend  themselves  to  analytical  treat¬ 
ment  or  easy  hand  computation.  It  is,  however,  feasible  to  calculate  them 
on  a  high  speed  computing  machine.  Thus  the  functions  F^,  F^,  . ...,  F^, 

0,  and  M  have  been  programmed  for  computation  by  the  EDVAC  in  the  Ballistic 
Research  Laboratories,  and  they  can  be  obtained  quickly  for  any  particular 
case  of  supersonic  flow  by  specifying  and  y  in  the  computing  machine  input. 

The  simple  form  of  the  expressions  in  Eqs.  (l.l)  suggests  the  practi¬ 
cability  of  a  compilation  of  the  coefficients  of  K  and  l/r  which  could  be 

s 

employed  as  a  labor  saving  means  of  determining  the  desired  quantities. 

To  this  end,  the  present  paper  provides  a  series  of  tables  of  shock  wave 
gradient  functions  versus  (3  for  y  =  1.4,  with  as  parameter. 

The  angle  £  is  of  special  interest  to  the  authors,  particularly  in 
calculating  a  flow  field  (velocity,  streamlines,  etc.)  from  given  inter¬ 
ferometric  density  data  (see  Ref.  5,  7>  and  8).  The  information  resulting 

from  the  measured  (3  and  K  reduces  the  uncertainty  in  density  near  a  shock 

s 

wave  which  is  inherent  in  the  interferometric  method. 

Sample  hand  computations  have  verified  that  the  formulas  used  in  this 
report  yield  the  same  values  as  formulas  given  in  the  references . 


9 


II.  DERIVATION  OF  GRADIENT  FUNCTIONS 


The  present  derivation  follows  the  one  indicated  in  Ref.  4.  In 
natural  coordinates  the  equations  of  isoenergetic  flow  of  an  ideal  non- 
viscous  gas  are 


/  -  /  \ 
u/p; 

op/ os  +  {±/ q;  oq/os  +  ow/on  +  e^sin  fcj/r  =  u 

(a) 

dp/ds 

+  pq  dq/ds  =  0 

ft) 

bp/bn 

o 

+  p  q  bQ/bs  =  0 

ft) 

h  +  q£ 

'  /2  =  constant 

(a) 

TdS  = 

dh  -  dp/ p 

(e) 

(2.1) 


where  e  =  0  for  two-dimensional,  =  1  for  axisymmetric  flow.  The  coordinates 


ovirl  v-i  o  O  vin  "1  nvi  rr4*  V-i  c 

CkXXlX  IX  CXJ.  ^  CU.  -l_^iXg>-UJU.C 


a  cma  r*  -4-  V-i  c*  4-  Yioowil  4  mci  ci  vi  a  "1  An  rr  4- Vnn  Viv\  rrAn  ol  fy«D  iorif 

OkX_WXXJ5  UIXV  .->1.1  iT-rmi  i  i  .  11-.  CXXXVA  r  i.  i  v  1 1  i  ^  UilVi  W  A  UUW  gWUCW. 


ory  of  the  streamline,  respectively;  q  is  the  speed  of  flow,  and  0  is  the 

angle  between  the  streamline  and  the  free  stream  direction  (x  axis);  p  is  the 

density,  and  p  the  pressure.  T  is  the  temperature;  M  is  the  Mach  number 

1/2 

(  =  q/a,  where  a  =  (dp/dp)  =  speed  of  sound). 


The  two  essential  equations  to  be  U6ed  here  are 

(l  -  M2)  (l/q)  dq/os  +  do/ai  =  -  c(sin  ®)/r  (a) 

(2.2) 

dq/dn  -  q  do/ds  =  -  (T/q)  dS/dn,  (b) 

obtained  from  Eqs.  (2.1)  and  the  definition  of  the  velocity  of  sound.  (The 
foregoing  equations  are  given  in  Ref.  9  and  10.)  Entropy  remains  constant 
along  streamlines. 


If  a  is  arc  length  along  the  shock  wave,  it  is  seen  from  the  relation 
b/ba  =  (b/bs)  ( bs/ba )  +  (d/dn)  (dn/da) 
that  at  the  shock  wave  ^2 

h =  ^  ^  =  Ks  h =  ^  cos  (p ' 9)  +  ^  sin  (p " e) 
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where  tan  p  and  K  are  the  slope  and  curvature,  respc  ctively,  of  the 
s 

shock  wave  (Fig.  l.l).  As  the  coordinates  are  defined,  curvatures  are 
positive  where  curves  are  concave  upward  and  negative  where  concave 
downward. 


By  means  of  Eq.  (2.3),  Eqs.  (2.2)  are  replaced  by  two  simultaneous 
linear  algebraic  equations  for  dq/ds  and  d0/8s.  The  coefficients  of  these 
equations  are  now  converted  into  the  desired  form;  namely,  that  of  Eqs. 
(l.l).  For  this  purpose  use  is  made  of  the  following  definitions  and 
relations  which  can  be  found  in  many  textbooks  and  reference  books  (the 
authors  employed  Ref.  6  mainly); 


M  > 

“1 


R 

K 


The  subscript  1  always  refers  to  free  stream  conditions, 


g  =  p^/'p  = 


(7  -  l)  m“  +  2 


*  =  Pj/P  =  (7  +  l)  /  jj 


q  =  a  q, 


where 


27m 


2  „  .  2  _,_2  „n1/2 


/  I  (7  +  l) 

L_ 

(7  -  1)] 


m 


A  =  (cos  0  +  g  sin  0) 

sin  (p  -  Q )  ~  (g  sin 

sin  9  =  Ul-  r)  sin 

IS 

t/^  =  g/i 
H/U±  =  A(Vg)1/2 


q.n 

s  =  (~) 

1 


( o  ),  N 


(2.5) 

(2.6) 

(2.7) 

(2.8) 


7(7  -  1)M, 


COS 

0  - 

A  \  / _ r>  \  f  a 

wy  =  p y/A 

(2.9) 

1 

J 

/A 

i  - 

(2.10) 

(2.11) 

(2.12) 

_ 

- 

(2.13) 

i°ge 

(j)  -  7  loge  (gOj 

+  S1 

After  differentiation  for  dq/dp,  d9/dp,  and  dS/dP  (noting  that 
6s/6s  =  0),  one  obtains  expression  for  F2,F^,  and  F^  given  by  the 


following  sequence  of  formulas; 
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L(m) 

=  g2  -  1  -  4g/(rm2  +  m2) 

(2.14) 

K(m) 

=  (1  - 

(2.15) 

■r  t  \ 

[g/(7  +  l)j  [jm2/  -  2(7  +  l)  +  2(7  - 

D/s] 

d  ^m; 

“  p 

f'y  _  1  ^  0  m  T  Y  m  ^ 

/  -  _  \ 
(,2.1b; 

fl 

=  g  [l(1  -  J)  -  K(1  -  gj] 

(2.17) 

f2 

-  -  K  f(L  -  g  +  l)/g  +  i2  i  L  (l  -  jjl 

(2.18) 

f5 

-  g  K(1  -  g) 

(2.19) 

=  K  (~(L  -  g  +  l)  +  L(l  -  j)] 

(2o20) 

F1 

=  [(cos  p)/A5]  (fx  sin2  p  +  f2  cos2  p) 

(2.21) 

F2 

r  -TI  0 

=  (cos  P)/A^J  (f^  sin^  p  cos  p) 

(2.22) 

F3 

=  [(sin  Pj/A2]  (f^  cos2  p  +  g  f3  sin2  p) 

(2.23) 

F4 

=  -  (g  f_  sin5  p  cos  P)/A2 

(2.24) 

bQ/be 

-  Fn  Kc  +  e  Fp  (l/r) 

(l/qj  bq/ba  =  F,  K_  +  €  Ft,  (l/r) 

•  ■  ~X  p  b  M- 

The  following  relations,  obtained,  from  the  flow  equations,  are  used 
to  find  £  and  i : 


P* 

P  =  P  7 

fit 


/  /  \  y 

«W' 


Vpt, )r  b 


/  /  \ 
^t/ 


(2.25) 


where  the  subscript  t  denotes  the  stagnation  value  (Pp^/Ps  =  Pp^/Ps  =  o), 


ana 


«v\  -  ‘v\  =  R  - i 


1/(7  -  i)  -7/(7  -  I) 


(See  Ref.  6) 
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-1/(7  -  1) 


-  [i  +  (r  -  1)  M^/a] 

*AX  *  [i  +  (r  -  i)  M^/a] 


-7/(7  -  l) 


By  means  of  Eqs.  (2.3)  ana.  (2.25),  dp/ds  and  dp/dn  can  be  expressed 
in  terms  of  dp/ds,  dp/dn,  and  other  quantities  already  determined  in  this 
report.  Substituting  these  expressions  into  Eqs.(2.l),  one  then  derives 
expressions  for  dp/ds  and  dp/dn;  and  putting  these  into  the  relation  for 
the  angle  between  the  density  contour  and  the  streamline,  namely. 


tan  £  =  -  (dp/ds)/ (dp/dn) 


(2.26) 


one  obtains 

tan 


.  _  ^  Ks  +  6  d/r) 

->  =:  T71 


Fc  K  +  e  F/ 
5  s  o 


TW  ' 


where 


F  = 

5 


1  COS  P 
7  +  1 


-  1 

gXm  _ 


AF, 


F6  =  -  af2 


(2.27) 

(2.28) 


A  positive  value  of  £  indicates  a  counterclockwise  rotation  from  the 
streamline  to  the  density  contour. 


The  angle  £.  between  the  Mach  number  contour  and  the  streamline  is 
determined  from 

tan  £  =  -  (dM/ds)/ (dM/dn)  =  -  (^l2/ds)/(dM2/dn),  (2.29) 

where  I  has  the  same  sign  convention  as  £.  Since 

M2  =  (pq2)/(7P), 

dM2ds  and  dM2 /dn  can  be  expressed  in  terms  of  the  derivatives  of  p  and  q, 
which  in  turn  are  expressible  in  terms  of  previously  determined  quantities. 
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Consequently  one  obtains 

K  +  e  F.  (l/'r) 
tan  £  =  r- 


frir^rT¥^n7?T 


,  where 


4  cos  p 

T 

a2  _ 

2g  (m2  -  l)2  sin2p  J> 

- - - - - p~~  > 

/  v  _l  1  ^  rr  m  flln  R 

n  — 

(y  +  l)  m2  1 

w  '  y 

(  0  ~zr\  ^ 

\  /  t  .*./  D  ™  r~ 

k 

J  (A5  7  l  m  /g)  +  2  A  sin  * 

+  ‘  .  2  . 

sin”  p 


s->  > 


Pl  + 


r, 


I 

l* 


A2  m2  i  cot  ft  I 

p  p 

g  sin  P 


2  +  (7  -  1) 


«2_2  , 
am  * 
o 

g  sin1-  p 


-1 


III.  DISCUSSION  CF  COMPUTATIONS 

Tabulations  of  the  functions  Fj,  F2,  ....,  F^,  9,  and  M,  calculated 
for  7  =  1.4,  are  presented  in  Section  4  for  values  of  ranging  from  1.1 
to  10.0*  The  functions  are  listed  in  order  of  decreasing  P,  from  p  =  90° 
to  a  value  somewhat  greater  than  p  =  sin  (l/M^) •  The  listed  values  of 
and  P  were  chosen  so  as  to  permit  the  determination  of  the  functions  to  a 
reasonable  degree  of  accuracy  for  any  and  p  by  means  of  graphical  inter¬ 
polation. 

Figure  3.1  shows  an  example  of  results  obtainable  from  these  calculations. 
The  right  hand  curve  shows  the  variation  in  slope  of  the  constant  density  lines 
along  the  shock  wave  of  a  supersonic  sphere  in  nitrogen  (shown  at  the  left) 
determined  from  shock  wave  coordinate  measurements  and  the  appropriate  shock 
gradient  functions.  For  reference  purposes,  the  left  hand  figure  shows  the 
sonic  point  and  "Crocco  point”  (defined  below),  and  indicates  the  location 
of  typical  values  of  rKR.  It  must  be  noted  that  the  accuracy  of  results 
like  those  in  Fisrure  3«1  Is  restricted  by  the  limited  accuracy  of  p  and  K 

W  S 
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DETACHED  SHOCK  WAVE. 


determined  from  experiment.  (See  Ref.  5  and  8  for  discussion  on  deter¬ 
mination  of  0  and  K  .) 

8 

A  point  often  referred  to  In  two-dimensional  flow  studies  (e.g. ,  Ref. 

3)  is  the  "Crocco  point”,  where  the  streamline  curvature  is  zero.  Figure 
3.2  shows  the  geometrical  conditions  which  must  exist  at  the  Crocco  point 
in  plane  and  axisymmetric  flows  for  7  =  1.4.  It  can  be  seen  that  in  two- 
dimensional  flow  the  Crocco  point  always  lies  below  the  sonic  point  on  a 
shock  wave  of  continuosuly  decreasing  slope.  The  two  points  approach 
coincidence  with  increasing  Mach  number.  On  the  other  hand,  in  axisymmetric 
flow  the  sonic  point  lies  below  the  Crocco  point  for  all  values  of  and 
rKg  shown.  By  continuity  there  would  be  values  of  these  two  parameters 
for  which  the  sonic  point  lies  above  the  Crocco  point. 

The  slope  of  the  sonic  line  is  a  quantity  of  considerable  interest.  It 
can  be  otained  from  the  tables  by  evaluating  the  (3  and  shock  gradient 
functions  corresponding  to  M  =  1  (e.g.,  by  graphical  interpolation);  it  is 
found  to  be  given  by 

(M1)  rKs  +  1 

1811  (9  +  ^  son.  -  Sp^TriT  +  B^TKJT’ 

where  B^,  B^,  and  B^  are  given  in  Table  3*1  for  7  =  1.4.  It  is  seen  that, 

formally,  rK  =  -  00  corresponds  to  two-dimensional  flow.  This  particular 
s 

case  is  plotted,  in  slightly  different  form,  in  Ref.  11. 

Figure  3.3  shows  plots  of  the  variation  of  ^son  >  the  angle  between 

the  streamline  and  Mach  number  contour  at  the  sonic  point.  For  two-dimensional 

flow  6  is  negative,  indicating  that  the  sonic  line  has  a  smaller  slope  than 

the  streamline.  For  axisymmetric  flow  the  figure  indicates  that  the  sonic 

line  can  have  a  greater  slope  than  the  streamline.  (The  streamline  curvature 

must  be  positive  in  this  case.  This  can  be  seen  from  Eqs.  (2.1b)  and  (2.1c) 

and  the  fact  that  dq/da  <  0  and  <3M/d a  y  0  along  a  shock  wave  with  continuously 

decreasing  slope.)  If  rK  =  -.2  at  the  sonic  point,  the  sonic  line  and 

s 

streamline  are  nearly  tangent  to  each  other  over  a  very  wide  range  of  Mach 
number . 
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Sonic  point 


Figure  3.2,  Conditions  for  zero  curvature  of  streamline 
id6/ds*Q)  at  shock  wave(Crocco  Point). 


i? 


(NOTE:  FOR  THE  rK$*0  CURVE,  0+£*£) 


FIGURE  3.3.  VARIATION  OF  ANGLE  BETWEEN  STREAMLINE 
AND  SONIC  LINE  AT  SHOCK  WAVE. 


The  values  of  M  sore  also  useful  in  finding  the  slopes  of  the 
characteristics  at  the  shock  wave  in  the  supersonic  portion  of  the 
flow  field. 


For  quick  reference,  Eqs.  (l.l)  are  restated  here. 
d©/ds  =F1Ks+  F2  (l/r) 

(l/qjK&j/ds)  =  F5  Ks  +  Fh  (l/r) 
tan  t  =  V1*.5  +  ri. 

fp;rrrr6 

tan  6  -  PJ<rKs>  +  F4 

W  *  p6 
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TABLE  5.1 


Sonic  Point  Table  (7  •  1.4) 


”1 

TJ 

^1 

R 

^2 

R 

12 

fl° 

r~ 

0° 

1.10 

10.05 

.793 

.301 

73.26 

1.40 

i  in 

x. 

6,15 

.175 

.464 

65.11 

6.31 

1.50 

4.78 

-  0.548 

.524 

62.26 

11.71 

1.75 

3.71 

-  1.40 

.546 

61.31 

17.77 

2.00 

2.78 

-  2.29 

.541 

61.50 

on  7A 
ZLZL  •  I  V/ 

2.50 

1.16 

-  3.91 

.517 

62.66 

29.65 

5.00 

-  0.51 

-  5.23 

.492 

63.77 

33.99 

3.50 

- 1.53 

-  6.27 

.474 

64.62 

36.82 

4.00 

r\  zz£ 

-  ^ 

7  nn 

-  f  .  \jy 

.458 

•a-w 

38.74 

4.50 

-  3.38 

-  7.70 

.451 

65.73 

40.06 

5.00 

-  4.04 

-  8.20 

.444 

66.08 

41.03 

6.00 

-  5.04 

-  8.92 

.455 

rzT  rr  Q 

UD.pU 

1.0  70 

7  nn 

|  •  v/v/ 

_  r  7n 

-  y  •  |  v 

-  9.38 

.428 

66.88 

43.16 

8.00 

-  6.19 

-  9.72 

.424 

67.09 

43.64 

9.00 

-  6.50 

-  9.92 

.420 

67.22 

44.03 

10,00 

Z'  rr/T 
-  Do  f D 

T  A  AA 
--LVJ  . 

.419 

£7 

U4*27 
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IV  -  TABLES  OF  SHOCK  GRADIENT  FUNCTIONS 
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